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We study the ground state properties of the Heisenberg spin-1/2 chain with ferromagnetic nearest- 
neighbor and antiferromagnetic next-nearest-neighbor interactions using two approximate methods. 
, One of them is the Jordan- Wigner mean-field theory and another approach is based on the trans- 

formation of spin operators to bose-ones and on the variational treatment of bosonic Hamiltonian. 
Both approaches give close results for the ground state energy and the magnetization curve at T — 0. 
£Nj , It is proved that quantum fluctuations change the classical critical exponents in the vicinity of the 

. transition point from the ferromagnetic to the singlet ground state. The magnetization processes 

' display a different behavior in the regions near and far from the transition point. The relation of 

the obtained results to the experimental magnetization curve in .R&2CU2M03O12 is discussed. 

m ; 

^ ■ I- INTRODUCTION 

U '. 

Lately, there has been considerable interest in low-dimensional quantum spin systems that exhibit frustration. The 
spin-1/2 chain with next-nearest-neighbor interactions (which is equivalent to a zig-zag ladder) is a typical model 
& ■ with the frustration. The Hamiltonian of the model is given by 

^ ■ N 

H = y^(JiS» ■ S n+ i + J2S W ■ S w+2 ) (1) 

Ch ■ n=l 

o : 

CJ , where S is the spin-1/2 operator, Ji and J2 denote the nearest (NN) and the next nearest neighbor (NNN) interactions. 

This model with both NN and NNN antiferromagnetic interactions (Ji, J 2 > 0) is well studied [J 0, 0, 0, 0, • 
\ There is a critical value J 2c = 0.2411Ji, which separates doubly degenerated dimer phase (at J2 > J 2c ) characterized 
■ by the excitation gap and the gapless spin- fluid phase (at J2 < J 2c)- Relatively less is known about the model 
with the ferromagnetic NN and the antiferromagnetic NNN interactions (Ji < 0, J2 > 0). Though the latter model 
has been a subject of many studies 0, 0, 0, EH the complete picture of the phases of this model remains unclear 
up to now. It is known that the ground state is ferromagnetic for J2 < ^p-. At J2 = ^p- the ferromagnetic state 
is degenerated with the singlet state. The wave function of the singlet state at J2 = ^p- is known exactly pdl IT^ . 
For J 2 > JAl the ground state is an incommensurate singlet. Controversial conclusions exist about the presence of a 

gap at J2 > -^r- ■ It has been long believed that the model is gapless 0, ^| but the one-loop renormalization group 
analysis indicates 0, Il4j that the gap is open due to a Lorentz symmetry breaking perturbation. However, existence 
of the gap has not been verified numerically 0] . On the basis of a field theory consideration it has been proposed E3 
1 ' that a very tiny but finite gap exists which can be not observed numerically. 

Besides a general interest in frustrated quantum spin models there is an additional motivation for the study of 
the model Q with J\ < 0. Really, the ferromagnetic NN interaction is expected to exist in compounds containing 
CuO chains with edge-sharing CuO^ units. The Cu — O — Cu angle in these compounds is close to 90°, and 
• • ■ usual antiferromagnetic NN exchange between Cu ion spins is suppressed. This means that the sign of J\ can be 
neg ative, while the NNN exchange is antiferromagnetic. Several compounds with edge-sharing chains are known 
[l6j . such as Li2Cu02, La§Ca%Cu2\0/t\, CCL2Y2CU5O1Q. However, in these compounds the antiferromagnetic long 
range order appears at low temperatures due to small interchain interactions. Recently, the crystal RU2CU2M03O12 
with edge-sharing chains has been synthesized and studied experimentally [TtI Il8j . Remarkably, no magnetic phase 
transition appears down to 2K, which testifies that interchain interaction is very small in this compound. Therefore, 
Rb-iC uoMo-tO 12 can be considered as an ideal model compound with the ferromagnetic (NN) interaction. According 
to [17[ it is described by the Hamiltonian {IJ with J\ ~ — 140-K" and ~ 0.4. 

One of the interesting peculiarities of the model is the dependence of the magnetization on the applied magnetic 
field at T — 0. The magnetization curve for Ji < is quite different from that for the case Ji > 0. It is characterized 



X 



'Electronic address: dmitriev@deom.chph.ras.ru 



2 



by a rapid increase (or even discontinuity) in the magnetization, if the external field exceeds a critical value. It is 
expected that the magnetization displays a true jump (the metamagnetic transition) when the NNN interaction J 2 is 

slightly larger than . This conclusion has been made on the basis of the exact diagonalization calculations for 

finite chains. However, the extrapolation of these results to the thermodynamic limit is rather difficult due to strong 
non-monotonic finite size effects. 

In this paper the ground state energy and the magnetization processes in the model Q with Jj < are studied 
using two variational approaches. One of these approaches is based on the Jordan- Wigncr transformation of spin- 1/2 
operators to the Fermi ones with the subsequent mean-field treatment of the Fermi-Hamiltonian. Another variational 
approach is applied to the Bose-Hamiltonian arising from a special transformation of spin operators to Bose ones. 

The paper is organized as follows. In Section II the details of the Jordan- Wigner mean-field approximation are 
presented. The zero-temperature magnetization process is studied and the form of the magnetization curve in different 
region of the parameter is found. For the case h — we have focused on the behavior of the model in the vicinity 
of the transition point from the ferromagnetic to the singlet ground state. The critical exponents characterizing 
this behavior are determined. In Section III we describe the variational method for the treatment of the considered 
Hamiltonian rewritten in bosonic form. Scaling estimates of the critical exponent near the transition point are 
presented in Section IV. In Section V we summarize our results. 



II. JORDAN- WIGNER MEAN-FIELD APPROACH 



It is convenient to represent the Hamiltonian of the spin- 1/2 chain with the ferromagnetic NN and the antiferro- 
magnetic NNN interactions in the form: 



H 



N 

£ 

n=l 



1 N 1 N 

(S„ • S„+i - -) + a^(S„ • S„+2 - j) -h^S* 



(2) 



71=1 



where a = pjj and h — 9 | J ^j B is the effective dimensionless magnetic field. The constant shifts in Eq.(J2J) secures the 
energy of the fully polarized state to be zero. 

Let us start from the classical picture of the ground state of Eq.J5J. In the classical approximation the spins are 
vectors which form the spiral structure in the XZ plane with a pitch angle tp between neighboring spins and all spin 
vectors are inclined towards the Y axis by an angle 0. The classical energy per site 



e(ip, 0) = i [1 — cos tp — a(l — cos(2<^)] cos 2 9 — ^ sin 9 



(3) 



is minimized by the angles 



4a 



(4) 



where 7 = a — j. 

The classical ground state energy is 



E c i _ I 2 
~N ~ ~2^ 7 



8 7 2 



(5) 



Following this picture we transform local axes on n-th site by a rotation about the Y axis by ipn and then by a 
rotation about the X axis by 9. The transformation to new spin-i operators r\ n has a form 



S„ = R y (tpn)R x (9)r} n 

where R y (ipn) and R x {9) are the operators of the corresponding rotations. 

Substituting (jfj) into J2J we obtain the transformed Hamiltonian in terms of the 77 operators 



H = H 1 +H 2 
H Y = Ne{cp,6)- 



H 3 

N 

£ 

n=l 



(6) 



(7) 



JlxVnVn+l + JlyVnVl+1 + JlziVnVn+l ~ 4) 
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N r ii^i 

J] J 2xVnVn+2 + J^X+2 + J2z(r)nVn+2 ~ ~ h sin 6 Yjtin ~ 7^ 



N 



H 2 = sin (9 X[sinv?«77f l+1 - r/X+i) - afa2<p(rfcr% + 2 - VlVn+2)] 



n=l 

N 



H 3 = ~ - cosip)(r) z n r} y n+1 + ^X+i) ~ a ^Y~0- ~ co ^W n vt+2 + ^X+2) 

n=l 



JV 



- sin if cos 0(VnVn+i - VnVn+i) ~ " sin 2^ COS @(?ln^ln+2 ~ VnVn+2)} ~ h cos 9 ^ rf n 



where 



J\ x = — cos ip, J\ y = — cos if sin 2 9 — cos 2 9, J\ z — — cos ip cos 2 — sin 2 9 

J2x = acos2tp, J 2y — a(cos 2</?sin 2 9 + cos 2 9), J 2z = a(cos ip cos 2 9 + sin 2 9) (8) 

The Hamiltonian (JJJ for h — (6* = 0) has been considered before in 0,0]. In Ref.fToj the ground state in the 
vicinity of the point a = j has been studied on the basis of the spin-wave theory taking (p by its classical value. It 
was shown that the transition at a = i from the ferromagnetic state to the singlet one is of the second order and the 
ground state energy is —Aj 2 for < 7 <C 1. In 0] the dependence of the ground state energy and the pitch angle <p 
on a has been found using the coupled cluster method. In particular, the ground state energy is proportional to 7 2 
for < 7 « 1 too. However, both cited approaches are not variational ones. Besides, the magnetization curve has 
not been studied. 

Our primary interest is how quantum effects alter the classical ground state structure. In this section we use the 
Jordan- Wigner transformation which converts the spin Hamiltonian Q into a model of interacting spinless fermions 



n — 1 ^ 

Vn = c n exp(i7r c+c 3 ) , jf n = - - c+c„ (9) 
The Hamiltonians H\ and H2 in Eq.(JZJ) are transformed to the Fermi-Hamiltonians having forms: 

N N N 

H lf = Ne((p, 9) + (hsin9 - J lz - J 2z ) ^ c+c„ + J lz ^ c+c„c+ +1 c„+i + J 2z ^ c+c„c+ +2 c n+2 

n—1 n—1 n—1 

1 N 1 N 

+ ^(Jlx + Jly) ^(C+C^+I + C+ +1 C„) + -{J 2x + J 2y ) ^(c+C„ +2 + C++2 C «)(! ~ 2c „+l c n+l) 
n—1 n—1 

1 N 1 N 

+ ^(Jlx ~ Jly) X]( C ™ C n+l + C «+1 C «) + l( J 2* - J 2!/) X]( C n C n+2 + C n+2 C n ){\ - 2c+ +l C n+1 ) (10) 
n—1 n—1 

AT 

H 2f = ^ sin6>^[sin99(c+c„ + i - c+ +1 c„) - asin2<p(c+c„ +2 - c+ +2 c„)(l - 2c+ +1 c„+i)] 

Z n=l 

We do not present here the Hamiltonian H3 f because it has very complicated form containing non-local interaction 
and, as it will be shown below, does not contribute to the energy in the mean-field approach. 

The next step of the approach is to treat the Hamiltonian Hf — H\f + H 2 f + H^f in the mean-field approximation 
(MFA). We use the approximation scheme proposed in pflj . The ground state wave function in this approximation 
has BCS-like form 



I*) = n(cos0 fe + sin^ C + C + fc )|O> 



fc>0 



where c£ is Fourier transform of the Fermi-operators c+. 

The expectation value of H 2 f over the function 'J is (H 2 f) = due to the fact that Sfe sm ^( c fe c fe) — an d 
^2k cos ^ ( c k c -k) = 0- The Hamiltonian H^f includes odd number of the Fermi-operators and, therefore, it is zero 



where 



in the MFA (i?3/) = as well. Thus, both Hamiltonians H 2 f and H^f do not contribute to the energy in this 
approximation. 

In the MFA four-fermions terms in Eas. are decoupled in all possible ways. After Fourier transformation the 
mean-field Hamiltonian takes a form 

Hmfa = NC + J2 a ( k )(4 c fc + c -k c -k) + J2 b ( k )( c t c -k + c -kCk) (11) 

fe>0 k>0 

C = e(^6)-(J lz + J 2z )n 2 + J lz (t 2 1 -g 2 1 ) + J 2z (t 2 2 -g 2 2 ) + (12) 
{■hx + J2 V ){nt 2 - 1[ - g\) + (J 2 x - J2 V )(ng 2 - 2tigi) 

a(k) = u + v\ cos k + v 2 cos 2k (13) 
b(k) = w\ sin k + w 2 sin 2k 

-{Jiz + J2z)(l - 2n) - (J 2x + J 2y )t 2 - (J 2x - J 2y )g 2 + hshx9 (14) 



u = 

V\ = 2^ lx ~ 2Ji z <i + 2(J 2x + J 2y )ti + 2{J 2x — J 2y )gi 

v 2 = {J 2x + J 2y )(^-n) -2J 2z t 2 

w\ = -{J\ x ~ Jiy) + 2{Ji z + J 2x + J 2y )g 1 + 2{J 2x - J 2y )tx 
w 2 = (J 2x - J 2y ){^ - n) + 2J 2z g 2 
The ground state energy, the one-particle excitation spectrum e(k) and the magnetization M = (S%) are: 

| = e(<p,6) + F( i p,9) 



M = I - -n) sind (15) 



where 



e(k) = ^Ja 2 (k)+b 2 {k) 
'1 

F(ip, 9) = {J lz + J 2z )n(l -n) + J lz (g 2 - t\) + J 2z (g 2 + t\ - 2t 2 ) + -(J lx + J ly )t x 

+ ^(Jix - Jx y )gi + (J2x + J2y){nt 2 - 1\ - g 2 ) + {J 2x - J 2y )(ng 2 - 2t 1 g 1 ) (16) 

Quantities n, ti, t 2l <?i, g 2 are the ground state expectation values, which are determined by the self-consistency 
equations: 

7T 

, , , f dk I a(k) 



(I 



t2 = (c+C„ +2 ) = - 

(i 

.91 = {c+c+ +l ) = - 





2tt V e(fc), 

dk a(k) cos k 
2^ e(fc) 

dk a(k) cos 2k 



2tt e{k) 

dk b(k) sink 
2^ e(k) 



, dk b(k) sin2fc 

32 = = - i T^HikT (17) 
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The solution of the self-consistency equations gives the minimum of the ground state energy in a class of a 'one- 
particle' wave functions at given angles ip and 9. We treat the angles ip and 9 as the variational parameters (not 
equal to their classical values). Thus, one should minimize the energy with respect to the angles ip and 6, solving the 
self-consistency equations for each value of ip and 9. This means that the proposed procedure remains variational one. 

To study the effect of the dimerization we added the staggered terms to the NN expectation values: 

(c+ Cn+1 ) = ii +*;(-!)" 

(c+c+ +1 ) = gi + g[(-i) n 

which leads to more complicated form of the Hamiltonian Hmfa and two additional self-consistency equations for t[ 
and g[. But the solution of the self-consistency equations with the minimization of the energy over ip and 9 gives 
t[ = g[ = 0. This means that the MFA does not indicate the dimerization in the system. 



A. Ground state energy near the transition point a — | 



4 



At first we consider the Hamiltonian at h = Q when 9 = and we are interested mainly in the behavior of 
the model in the vicinity of the point a = 4, where the transition from the ferromagnetic to the singlet ground state 
occurs. At h = the ground state energy per site is: 

|=e(^,0) + F(^0) (18) 

where F(ip, 0) is a quantum correction to the classical part of the energy. 

The analysis of the solution of Eas. (|17(l shows that in the case < 7 ^ 1 and ip <C 1 the functions e(ip, 0) and 
F(ip, 0) have forms: 



2 4 

2 32 

„4 



= -^+ A ^ ( 19 ) 
The coefficient A w 0.0195 is determined by the numerical solution of Eas. (|17fl . 

Substituting (|19fl into i|18|) and minimizing E with respect to ip we obtain the leading terms in 7 for the angle ^(7) 
and the energy E(j): 

ip = 2.33I7T1 
E 12 

- = -1.585 7 - (20) 

The numerical solution of self-consistency equations l|17fl confirms this dependence E(j) at < 7 <C 1 (see Fig^). 

Since the MFA is the variational approach, the found critical exponent for the ground state energy (3 — ^ gives upper 
bound for the exact critical exponent (3 < y- and certainly it is less than f3 = 2 obtained in classical approximation. 
We note that the use of more elaborate methods 0, 0] does not change the exponent j3 = 2 and change the numerical 
prcfactor a in the ground state energy E = —aj 2 only. This means that at present the MFA gives the best estimate 
of the ground state energy in the region < 7 <C 1- 

In Ref.[22] the model (j2J (at h = 0) has been studied using the Jordan- Wigner transformation and the mean-field 
theory. However, the pair correlations of type (c+c+ } have been neglected. It was conjectured in j2^ that the ground 
state at \ < a < 0.25854 is not singlet but has a nonzero total spin. We have found no evidence of this fact and the 
numerical diagonalization of finite chains does not confirm this prediction as well. 

As was written above, in our mean-field treatment the pitch angle ip is considered as the variational parameter. 
It should be noted that in the MFA a helical (spiral) long range order (LRO) exists as well as in the classical 
approximation. The helical LRO is characterized by the angle ip, though its value differs from the classical one. Of 
course, this LRO is destroyed by those quantum fluctuations which are ignored in the MFA. At the same time the 
obtained helical LRO is regarded to an incommensurate behavior of the correlation function and ip can be identified 
with the momentum g max at which the static structure factor takes a maximal value. 
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FIG. 1: The dependence of the ground state energy of the model ^ on a in the MFA (short-dashed line) and the boson 
variational approach (long-dashed line). Solid line corresponds to the energy given by Eg. 1201 . 



B. Magnetization curve 

In this section we consider the magnetization processes in the model (0- The behavior of the magnetization in 
the region of the magnetic held close to saturation is of a particular interest. We note that the determination of 
the saturation field h s (a) at which the transition to saturation takes place is generally not a simple problem. In the 
model with both antiferromagnetic NN and NNN interactions this held is equal to the energy of an one-magnon 
state on the ferromagnetic background (all spins up), i.e. h s = —Ex. It is not the case for the model ©■ The specific 
property of this model is the presence of bound magnon states (complexes of two or more flipped spins) arising as 
a result of magnon-magnon attraction. When the held h is reduced below the saturated value the number of spins 
flipping simultaneously is more than one. Then the saturated field is determined by the condition 

/j s = max|^| (21) 

where E n is the energy of n-magnon state. 

There are two possible scenarios of the behavior of the magnetization curve close to saturation. In the first case 
Eq.l|2ip is satisfied at n* ~ o(N) and the magnetization is a continuous function at h — > h s . In the second case 
n* ~ N and the magnetization jumps at h = h s from M* — {o ~ jr) to M = ^, i.e. the metamagnetic transition 
occurs. It means that the ground state energy as a function of magnetization e(M) has a negative curvature and M* 
is determined by the Maxwell construction. We note that the numerical estimate of h s and n* especially at < 7 <C 1 
is very difficult because the number of spins flipping simultaneously at M — > k increases when a — > 4 0, . 

We calculate the magnetization for the model J5J using the MFA at h > 0. In this case there are two variational 
parameters, ip and 6. We show the ground sate energy per site e as a function of h for two values of a in Fig|3 
together with the energy of the ferromagnetic state = —\- The saturation held h s is determined by the crossing 
between the variational energy e(h) and &f ■ The dependence e[K) shown in Fig [3 demonstrates two different types 
of the behavior. For a < 0.38 (Fig[2Ji) e(h) and ep(h) have different slopes at the crossing point. At h = h s the 
magnetization jumps from M* — §^\ h > to ^ as it is shown in FigEt- It means that maximum in Ea. (|21|l is reached 
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FIG. 2: The dependence of the ground state energy of the model (2) on magnetic field e(h) for (a) a = 0.3 and (b) a = 0.4. 




for macroscopic n and the metamagnetic transition takes place. The dependence of the saturated field on a is shown 
in Fig0] This dependence is consistent with those obtained in [l9l | by using numerical diagonalizations of finite chains. 
The magnetization curve at < 7 <C 1 can be found analytically from the self-consistency equations l|I7(l . The energy 
per site e(M) at M -C 1/2 has a form 

e(M) = -1.5857^ + 3.69 7 ^Af 2 + 0( 7 ^) (22) 

Performing the Maxwell construction we obtain 

M* = 0.4297? 

h s = -3.177^ (23) 

As follows from Eq. H21|l h s in Eq. (|23|l is the binding energy per magnon of multimagnon complex, ft is interesting 
to compare this value with the binding energy of two-magnon complex which is — 727 3 at < 7 <C 1 ED- It is 
evident that interaction of a macroscopic number of magnons essentially lowers the binding energy. 



8 




FIG. 4: The phase diagram of the model @ in the MFA in the (h, a) plane. 



The magnetization M(h) at < 7 <C 1 is 



M(h) = 0.1367 7 /&, h < h s 
1 

2' 



M(h) = i, /i>/i s (24) 



According to Eq.J23J) the susceptibility x is 

X = 0.1367 - ^ (25) 

As follows from Eq.()24p and as shown on Fig[3]the magnetization is a linear function of h at h — > 0. This testifies 
the absence of the gap in the spectrum. As it is noted before, the subtle question about the gap in this model is still 
open. Numerical calculations show 1 /N behavior for a gap |2j , while the one- loop renormalization group indicates so 
tiny exponentially small gap [TBi] that it can not be observed numerically. The proposed version of the MFA does not 
predict the dimerization and the gap for the model @. However, we are not sure that the mean-field approach is an 
adequate tool to answer so subtle question. 

For a > 0.38 the curve e(h) (FigEt>) has a cusp at some critical magnetic field h c as a result of the jump of 
variational parameters minimizing energy from ip, 9 < ^ to tp = 9 = 5. At h — h c the magnetization jumps from 
Mi = _ 9^0 tQ M ^ = _ d^h± _ According to Eq.@ e(h) at h > h c ((p = 9 = f ) is the ground state 

h—h c — h—h c + 

energy of the Hamiltonian 

N N N 

H=- " 4) + <* E( S « S «+2 - 4) - k E S n (26) 

n— 1 n— 1 n— 1 

The curve e(/i) in FigEt is tangent to ep(h) at the crossing point. The saturation field h s and the magnetization 
M (h) near h s are determined from the self-consistency equations l|17|l with ip = 9 = 5 
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M ~ 2 (l + 2a)(4 a 4 ( +15a3 + lla2-l) (/ls k) (27) 

where 

4a 2 + 2a - 1 

E 2 = — 28 

1 + a 

It is interesting to note that E 2 given by this equation is the exact energy of the bound two-magnon state with a 
momentum q = ir of the Hamiltonian i|26[l . The wave function of this state has a form 



TV 

with a correlation length 



J2(-l) n+l e- m S-S- +2l+1 |TTT - T) (29) 



r i=lni±^ (30) 



a 

Remarkable, that this function is the exact one of the two-magnon state of the Hamiltonian J5J as well because 



E 



(S*SZ +1 + S%S v n+1 ) M = (31) 



The state Ij29(l has the lowest energy in two-magnon sector of the Hamiltonian J2Jl at a > 0.38 0- 

Thus, the maximum in Ea. (|21|l is reached at n — 2 if e(h) has the form similar to that shown in Fig^- This fact 
confirms the observation based on finite-size results that in the region of a rather far from a = j only the flipping 
of pairs of spins participate in the magnetization process near the saturation. 

The dependence M(h) for two typical cases is shown in FigOU At a > 0.38 the MFA gives the saturation field exactly, 
though the linear behavior of M(h) in Ea. (|27|l is in a contrast with the expected square-root critical dependence. 

It was determined above that the presence of the cusp in the dependence e(h) is the result of the jump of variational 
parameters. At h < h c the angle (p is incommensurate while at h > h c it corresponds to the commensurate phase. 
Therefore, the magnetic field h c can be identified with the critical field at which the incommensurate-commensurate 
transition takes place. The phase diagram of the model (|2J) in the MFA is shown in FigQ] The commensurate phase 
lies between the incommensurate and the ferromagnetic phases. 

In the limit a — > oo the Hamiltonian J2J describes two independent Heisenberg chains. The numerical solution of 
the MFA equations indicates that in this limit the width of the commensurate phase shrinks as I /a. In the limit 
a — > oo the pitch angle <p tends to ir/2 for all values of h, while the canted angle 6 changes from = at h = 
to 9 = it/2 at h = h s . This means that the MFA correctly describes the limit a — > oo showing the incommensurate 
phase at < h < h s . 

The MFA shows the jump in magnetization at some critical field h c . We are not sure whether a true magnetization 
curve has (or not) the jump. However, we believe that the jump obtained in our approximation testifies the singular 
behavior of the magnetization at some critical field. A plausible reason of such singularity consists in the following 
physical picture. As was shown above at a > 0.38 two magnons attract one another and form the bound two-magnon 
state (|29H with the correlation length £ (|30|) . When the magnetic field is slightly lower than the saturated field h s 
and a number of magnons is small, the ground state can be represented with high accuracy as a product of the bound 
magnon pairs, weakly interacted with each other. The bound pairs of magnons start to feel each other when the 
mean distance between bound pairs becomes of the order of the correlation length of the pair £. This happens at 
some critical magnon density n c ~ For larger density n > n c one can not consider magnon pairs as independent 
quasi-particles and the above picture is destroyed. Certainly, these speculations can be well justified in the case a 3> 1 
only, when the critical density is small n c ~ 1/a. 

Possible singularity in the magnetization is related to the magnetization curve in RbiCuiMo-iOn . The experimental 
magnetization at T = 2.6K demonstrates a sufficiently sharp increase up to M ~ 0.4 at B ~ 14 Tesla followed by 
a gradual increase to the saturation. In the light of above mentioned one can assume that the sharp change in the 
behavior of M(h) at B ~ 14T is connected with the transition from the incommensurate regime to the commensurate 
one. Taking the values of J\ and a estimated for Rb 2 C '112M 03O '12 (a — 0.4, J\ ~ — 140JT) (13 we obtain B c = 13. 6T. 



III. BOSON VARIATIONAL APPROACH 



Another technique which can be used for an approximate analysis of the model (J2Jl is the boson variational approach 
(BVA) proposed in Ref. 23]. This approach is based on Agranovich-Toshich boson representation of spin s = 1/2 
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operators 24]: 

s,, 



2 



= bn - ■ . (32) 
2V^&n 

where 6+ are the Bose operators. This transformation preserves all commutation relations for the spin operators 
and does not involve 'nonphysical states', because all states with odd number of bosons correspond to the same spin 
state II) and the states with even number of bosons correspond to the spin state |T). The latter implies that the 
transformation i|32|l is not unitary and each energy level of spin model corresponds to infinite degenerated level of 
boson Hamiltonian. However, this transformation allows to estimate the ground state energy of the spin model. 

The Hamiltonian of the considered spin model rewritten in the bosonic form is treated by the variational function 
in the form pf: 

|*)=cxp l^A(n)b+bf +n ) |0 b ) (33) 



1. 1 1 



where |0;,) is the boson vacuum state corresponding to the ferromagnetic spin state with all spins pointing up and 
the function A(n) is chosen by the condition of minimum of the total energy. Therefore, in contrast to spin-wave 
theory, this approach is variational. The procedure of calculation of the variational energy with the wave function 
(|33JI and energy minimization over function A(n) was developed in 23], where the approach was successfully applied 
to construction of the ground state phase diagram of the frustrated 2D Hciscnberg model. 

We have applied the above approach to the rotated spin Hamiltonian Q and the rotation angles ip and 9 was used as 
variational parameters, as was in the MFA. The contribution to the energy of the parts H2 and H3 of the Hamiltonian 
in this approach is also zero as in the MFA. We do not present here cumbersome details of calculations, because 
they are identical to those in Ref.|2^] and give only the final results. 

The approach shows very similar behavior of the ground state energy as in the Jordan- Wigner MFA. That is, for 
small < 7 <C 1 both approaches give the same critical exponent for the ground state energy and for the parameter 
if. Moreover, the numerical estimates of the ground state energy in both approaches are also very close (see Fig^). 
The behavior of the magnetization curve in the BVA are also very similar to that in the MFA (see FigEJt - The perfect 
coincidence of physical properties of the model predicted by the MFA and the BVA looks somewhat surprising, because 
of different nature of these approaches. 

We note, that the potential advantage of the boson variational approach consists in its applicability to any dimension 
lattices and possible modification to any spin value S. 



IV. SCALING ESTIMATE OF THE CRITICAL EXPONENT NEAR THE TRANSITION POINT a = 1/4 

As was shown in previous sections the MFA and the BVA gives the ground state energy E^) ~ 7 12 / 7 near the 
transition point a — 1/4. Since both approaches used are variational, the value (3 = is upper bound for the critical 
exponent (3. This implies an important and strict fact that quantum fluctuations definitely change the classical critical 
exponent f3 = 2. 

Nevertheless, there are some reasons to believe that true value of the critical exponent is (3 = |. The arguments in 
favor of this conjecture are based on the following consideration with use of the results of numerical calculations of 
finite-size chains 0, Hi], |25j . 

The ground state of the model J2Jl is the singlet at a > \ for any even TV. For cyclic chain there are number of level 
crossings of two lowest singlets with momenta q = and q = it. For example, for N = Ak there are k — 1 crossing 
points j < ai(N) < a2(N) < . . . < ak-i(N) 0,I3- This fact means that in the thermodynamic limit the ground 
state is at least two-fold degenerated at a > j. Let us determine the dependence of a± (actually, 71 = a\ — j) on 
N . It is known this dependence defines the scaling parameter in one-dimensional models, where the crossings of two 
ground state levels occur. At a — j (7 — 0) the ground state in the sector 5 = has momentum q = 7r, while the 
first excited state has q = and the excitation energy AE ~ N~ 4 25]. The first order correction to the energy of 



11 



these states in 7 is 

^1(3) (7) =7E(*1P) l S « ' S «+2l *1(2)> (34) 
n 

where ^1 and \&2 are the wave functions of the ground state and the first excited singlet at a = j. These two states 
have a spiral ordering at N — > 00 with a period of spirals N and ^ respectively 0, ■ According to this fact the 
two-spin correlation functions at N — > 00 are 

<*1 |S„ • Sn+il^x) = | cos ^ 
1 47r/ 

<* 2 |S„-S, i+/ |vI/ 2 ) = jcos — (35) 

The accuracy of the above equations for I = 2 is of the order of 0(N~ 3 ) flU |25|. 
The value 71 is determined from the condition 

SExM - 5E 2 (ji) = AE (36) 

which gives 

7l ~ N- 3 (37) 

Therefore, the scaling parameter of the model in the vicinity of the transition point a = i is x = jN 3 . It means 
that the perturbation theory in 7 contains infrared divergencies and the ground state energy has a form 

E(7) = lf(x) (38) 

where the scaling function f(x) at x — > is given by the first order correction (|34() . In the thermodynamic limit 
(x — > 00) the behavior of f{x) is generally unknown, but the condition E ~ iV at iV — » 00 requires 

£?( 7 ) ~ N'f (39) 

According to Eas. l|38|) - (|39|l (3 — |. This value is close to /3 = which indicates high quality of the variational 
approaches used. A possible reason of a discrepancy between the variational and scaling estimates of the exponent 
may be related to the fact that in the variational approaches the terms H2 and H3 in Eq.J7I) are irrelevant. It can 
be also expected that the true dependences in Eas.l|2U |) .l|23 |l at small < 7 <C 1 are 

y.~ 7 3, M*~ 7 i, h s ~ji, x~7"3 (40) 



V. SUMMARY 



We have studied the frustrated spin-^ Heiscnberg chain with the NN ferromagnetic and the NNN antiferromagnetic 
exchange interactions using two different variational approximations: the MFA and the BVA. The first step of both 
approaches consists in rotation of the coordinate system on the pitch angle ip and the canted angle 9, which are not 
equal to their classical values and used as variational parameters of the approaches. Then, in the MFA the rotated spin 
Hamiltonian is mapped into the model of interacting spinless fermions by means of the Jordan- Wigner transformation. 
The latter is treated in the mean-field approximation with the inclusion of superconductor-like correlations. Within 
the BVA we use the Agranovich-Toshich boson transformation of spin operators and the variational treatment of the 
bosonic Hamiltonian. A remarkable fact, despite the difference of these approaches they give quantitatively close 
results. 

The variational approaches used allowed us to estimate the critical exponent of the ground state energy (3 in the 
vicinity of the transition point from the ferromagnetic state to the singlet one. Both approaches give j3 = which 
differs from the classical value (3 = 2. Since approaches used are variational, we have established an important and 
strict fact that quantum fluctuations definitely change the classical critical exponent. Using the results of finite-size 
calculations we presented also some scaling arguments in favor of the critical exponent j3 = | . This value is close to 
(3 = y- which indicates high quality of the variational approaches used. 

The behavior of the magnetization curve is different in parameter regions a ^ 0.38 and a ^ 0.38. In the region 
a ;$ 0.38 the metamagnetic transition to saturation takes place. At a > 0.38 the magnetization increases with h until 
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the field reaches the critical value h c , where the magnetization jumps from Mi to M.% < \. This magnetization jump 
accompanies the jump of the pitch angle <p and canted angle 9. At h < h c both angles are incommensurate, while 
at h > h c they correspond to the commensurate phase ip = 6 = w/2. Therefore, we associate this jump with the 
incommensurate-commensurate transition induced by the magnetic field. 

We believe that our approaches correctly predict the existence of the incommensurate-commensurate transition at 
some critical field h c and a > 0.38. This transition must accompany some singular behavior of the magnetization 
curve, though we are not sure whether the true magnetization has the jump at h = h c . 

Both approximations used yield the exact value of the saturation field at a > 0.38, though the field dependence near 
the saturation does not correspond to the expected universal square-root behavior. We note that the magnetization 
behavior near the saturation in the region of sufficiently large value of a is described by the model J2HJ), which looks 
simpler than the model J5J). In this respect an accurate study of the model H26[) is of a particular interest. 

We note that the magnetization curve in Fig.® resembles (apart from the magnetization jump) that observed in 
Rb2Cu2Mo^O\2- This resemblance allows us to provide a qualitative explanation of the peculiarity of the magne- 
tization process in this compound. The experimental curve is characterized by a sharp change of the susceptibility 
at the magnetic field B ~ 14 Tesla. We assume that this peculiarity in the magnetization behavior is related to the 
crossover between the incommensurate and the commensurate states. 

This work was supported under RFBR Grant No 03-03-32141. 
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